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Chung  In  Um 

Department  of  Physics,  College  of  Science 
Korea  University,  Seoul  136-701,  Korea 

The  Schrodinger  equation  is  used  to  exactly  evaluate  the  propagator,  wave 
function,  energy  expectation  values,  uncertainty  values  and  coherent  state  for  a 
harmonic  oscillator  with  a  time -dependent  frequency  and  an  external  driving  time- 
dependent  force.  These  quantities  represent  the  solution  of  the  classical  equation 
of  motion  for  the  time -dependent  harmonic  oscillator. 

I .  Introduction 

It  is  well  known  that  an  exact  solution  of  the  Schrodinger  equation  is 
possible  only  for  special  cases.  For  this  reason,  approximate  methods  are  needed. 
Exact  solutions  provide  important  tests  for  these  approximate  methods  and  for 
various  models  of  physical  phenomena.  In  general,  the  solution  of  the  Schrddinger 
equation  for  explicit  time -dependent  systems  has  met  with  limited  success  because 
of  analytical  difficulties,  although  progress  has  been  made  during  the  past  three 
decades.'*'  **  Camiz  et  al^  have  obtained  the  wave  functions  of  a  t  ime  -  dependent 
harmonic  oscillator  perturbed  by  an  inverse  quadratic  potential,  using  the 
Schrodinger  formalism  and  a  generating  function.  Further,  Khandekar  and  Lawande^ 
have  evaluated  the  exact  propagator  and  wave  function  for  a  time -dependent  harmonic 
oscillator,  both  with  and  without  an  inverse  quadratic  potential,  using  Feynman 

g 

path  integrals.  In  addition,  Jannussis  et  al  have  calculated  the  propagator  for 

several  quantum  mechanical  systems  with  friction. 

9 

In  a  previous  paper,  we  have  evaluated  the  propagator,  wave  function,  energy 
expectation  values,  uncertainty  values  and  transition  amplitudes  for  a  quantum 
damped  driven  harmonic  oscillator  by  using  path  integral  methods.  Also,  we  have 
obtained  the  coherent  state  for  the  damped  harmonic  oscillator*"^*  and  calculated  the 
propagator  for  coupled  driven  harmonic  oscillators.*'*' 
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In  this  paper  we  discuss  the  exact  quantum  theory  of  a  forced  harmonic 
oscillator  with  a  time -dependent  frequency.  In  Sec.  II  we  evaluate  the  propagator 
using  the  Schrodinger  equation  and  path  integral  methods,  and  in  Sec.  Ill  we 
calculate  the  wave  functions  using  the  propagator.  In  Sec.  IV  we  define  the  energy 
operator  and  calculate  energy  expectation  values.  In  Sec.  V  we  obtain  the 
uncertainty  values.  In  Sec.  IV  we  determine  the  coherent  state  and  its  properties. 
Finally,  in  Sec.  VII  we  present  results  and  a  discussion. 

II.  Propagator 

We  consider  a  system  whose  classical  Hamiltonian  is  of  the  form 

H  “  2M  p2  +  2  w2(t)  ^  '  f(C)x  1  (2.1) 

where  x  is  a  canonical  coordinate,  p  is  its  conjugate  momentum,  u(t)  is  a  frequency 
as  a  function  of  time,  M  is  a  positive  real  mass,  and  f(t)  is  an  external  driving 
force.  The  Lagrangian  corresponding  to  the  Hamiltonian  (2.1)  is 

L  -  ^  Mx2  -  2  Mu2(t)x2  +  f ( t)x  .  (2.2) 

Here,  the  Hamiltonian  H  and  Lagrangian  L  depend  on  time.  The  classical  equation  of 
motion  for  our  system  is 


d"  2.  ,  1  . 

— “  x  +  a>  (t)x  -  -  f(t) 

dt“ 


(2.3) 


For  the  case  where  w(t)  -  (constant),  the  solution  of  Eq .  (2.3)  represents 

harmonic  motion;  otherwise,  it  is  difficult  to  evaluate  the  exact  solution. 

The  path  integral  formulation  of  Feynman  provides  an  alternate  approach  to 
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solving  dynamical  problems  in  quantum  mechanics.  In  this  approach,  the  usual 
Schrodinger  equation  is  replaced  by  the  integral  equation 

i/>(x,t)  -  J  dx’  K(x,t;x't')  0(x',t')  (t  >  t')  (2.4) 

with  the  initial  condition  0(x,t)  -  ^>(x',t).  Here,  ^>(x,t)  is  a  wave  function  and 

K(x,t;  x',t')  is  a  propagator.  The  propagator  K(x,t;  x',t')  is  defined  by  the  path 

,12 

integral 

f(x,t)  N-l 

K(x,t;  x',t')  -  iim  II  dx.  exp[n  S(x,t;  x',t')] 

J(x',t')  j-1  J  ^ 


(2.5) 
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where  che  integration  is  over  all  possible  paths  from  the  point  (x',t')  to  the 
point  (x,t),  and  S(x,t;  x',t')  is  the  action  defined  as 


S (x , t ;  x',t') 


dr  L(x,x,r) 


(2.6) 


For  a  short  time  interval  «,  substitution  of  Eqs .  (2.2)  and  (2.6)  into  Eq.  (2.5) 
gives  the  normalizing  factor  and  the  usual  Schrodinger  equation: 


A. 

J 


(2.7) 


yj>  +  J  Mw2(t)  x2  rj>  -  f(t)  x  yj> 

2M  3x“ 


(2.8) 


Since  K(x,t;  x't')  can  be  thought  of  as  a  function  of  the  variables  (x,t)  or  of 
<x',t'),  it  is  a  special  wave  function,  and  it  satisfies  Eq.  (2.8): 

2  „2 

i#  K(x,  t;  x',t')  -  -  ^  2—^  K(x,t;x',t')  +  ^  Mw2(t)  x2  K(x,t;  x'  ,  t') 

3x 


K(x , t ; x' 


-  f(t)  x  K(x, t;x' , t' ) ,  (t  >  t') 

2  „2 

t')  ~  '  m  2 — 2  K(x-t:  x,t'>  +  2  x'2 

1  5x'~ 

x  K(x,t;x',t')  -  f(t')  x'  K(x,t;x',t') 


(2.9) 


(f 

Because  the  Lagrangian  is  quadratic,  the  propagator  has  the 
K(x,t;  x't')  -  exp(a(t, t' )x2  +  b(t,t')xx'  +  c(t,t')x'2 


>  c) 


form 
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+  g(t,t')x 


(2.10) 


+  h(t, t' )x'  +  d(t, t'  )  ] 


(2.11) 


where  from  Eqs.  (2.9)  and  (2.10)  we  can  easily  deduce  that  the  coefficient  of  the 
third  and  higher  powers  in  x  is  zero. 

Substituting  Eq .  (2.11)  into  Eqs.  (2.9)  and  (2.10),  we  obtain  the  differential 
equations 
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d_  21H  2  M  2 .  . 

51  a  “  M  3  +  2iH  “  (t) 


(2.12) 


^  (bx'+g)  -  ^  a(bx'+g)  +  J  f(t) 


(2.13) 


^  (cx'2+hx'+d)  -  &  a  +  ^  (bx'+g)2 


(2.14) 


3F  c  •  5?  c  +  2*  “(c  > 


(2.15) 


^7  (bx+h)  -  g  c  (bx+h)  +  “p 


(2.16) 


d  ,  2  _M_  ,,  ,  .  2  &_ 

—  (ax  +gx+d)  -  Jg  (bx+h)  +  gj  c 


(2.17) 


Equations  (2.12)  and  (2.15)  are  nonlinear  equations.  For  the  case  where  w(t)  -  wq, 
a  solution  is  easily  found,  but  in  other  cases  it  is  difficult  to  find  an  exact 
solution.  If  q(t)  obeys  the  differential  equation 


- -  q(t)  +  (j  (t)  q(  t)  -  0 

dt 


(2.18) 


then  the  solutions  of  Eqs .  (2 . 12) - (2 . 14)  are 


,  .  _  IM  al£l 

1  '  q ( t) 


(2.19) 


b(t)x'  +  g(t) 


i  1  rf‘ 

K  q  -J 


ds  f(s)q(s)  +  b 


(2.20) 


2 

b  rt 


:(t)x'2+  h(t)x'  +  d(t)  -  lnq'h  + 

2LlW  J  q(s)2  lM*  J  q ( s ) 


x  dp  f (p)q(p)  + 


±_  rc  ds  rs 

L>IM  J  ,  ,2  J 
J  q(s)  J 


dp  f (p)q(p)  dr  f(r)q(r)  +  d 


(2.21) 


where  b  and  d  are  constants  of  integration  and  do  not  depend  on  t,  and  the 
o  o 

solutions  of  Eqs.  (2 . 15) - (2 . 17)  are 


c(t')  - 


_M_  flitll 
2i#  q(t') 


(2.22) 


b ( t ' ) x  +  h( t' ) 


-ife'f 


ds  f(s)q(s)  +  b' 


(2.23) 
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a(t')x2  +  g(t')x  +  d(t')  -  Inq(t')  ^  + 


ib'  Pt' 


ib'  r  t 


•LU  ,  iU  , 

o  ds  _ o  ds 

2#M  J  .  .2  MH  J  .  . 

J  q(s)  J  q(s) 


f 


x  dp  f (p)q(p)  + 


—  f  —  f 

2iM»  J  q(s)2  J 


dp  f(p)q(p)  J  dr  f(r)q(r)  +  d^  , 


(2.24) 


where  b^  and  d^  are  constants  of  integration  and  independent  of  t'  .  Since  only  t 
is  a  variable  in  Eqs  .  (2 . 19)  -  (2 . 21)  ,  we  have  suppressed  t'  in  a(t,t'),  b(t,t'), 
etc.,  and  we  have  similarly  suppressed  t  in  Eqs.  (2 . 22) - (2 . 24) . 

In  polar  form  we  may  write 

q(t)  -  f7<t)e17<t)  ,  (2.25) 


where  r;(t)  and  7(t)  are  real  quantities. 
r/(t)  -  T7(t)  72(t)  +  cj2(t)  r?(t)  -  0 
2r/(t)y(t)  +  rj(t)  y(t)  -  0 

n2(t)7(t)  -  n 


From  Eqs.  (2.18)  and  (2.25)  we  note  that 

(2.26) 

(2.27) 

(2.27' ) 


•..here  the  constant  Q  is  a  time-variant  quantity.  From  Eqs.  (2.26)  and  (2.27),  we 
find  another  form  for  the  solution  of  Eq .  (2.18)  as 


q(t)  -  rj(t)  sin(7-7'  ) 


(2.28) 


q(t'  )  -  ri(t')  sin(7-7'  ) 


(2.29) 


where  7  -  7(t)  and  7'  -  7(t'). 

Substitution  of  Eq.  (2.28)  into  Eqs.  (2.19)  and  (2.21)  gives 


a(t)  -  ^  [J  +  7  cot(7-7') 


ib 

b(t)x'  +  g(t)  -  - — 7 - 7  +  - — 

H»7sin(7-7)  K»?s 


-1 - f‘ 

in(7-7')  J 


ds  r?(s) f(s)  sin[7(s)-7'] 


(2.30) 

(2.31) 


•x2 

2  -W  .u  lb0 

c(t)x'  +  h(t)x'  +  d(t)  -  Ln[r)  '  sin  '(7-7')!  +  - 


2hftM 


co t ( 7 -  7 ' ) 


6 


- -  fC 

n(7*7')  J 


i#flM  sin(7 

1 


ds  f(s)^(s)  sin(7(s) -7' ) ] 


+  2iJ<flM  sin(7 


^7y  ds  |  dp  f  (s)f  (p)r?(s)r?(p) 


x  sin(7(s)-7' ]sin(7(p)-7') 


(2.32) 


Furthermore,  substitution  of  Eq.  (2.29)  into  Eqs .  (2.22)  and  (2.24)  gives 


c(t')  “  ™  t-  +  7*  cot(7-7')] 

2J»  f? 


(2.33) 


b(t' )x  +  h(t* ) 


+  .  1 


_L 


i{<ij '  sin(7*7 '  )  i)tff?sin(7 


-f' 


ds  f  (s)f7(s)sin[7-7(s)  ] 


(2.34) 


9  _u  .k  ibo 

a(t')x  +  g(t')x  +  d(t')  -  ln[rj'  sin  (7-7')]  +  0WnM  cot(7*7') 


2hOM 


ib 

hflMsin(7 


rC* 


LU  pi. 

in(7-7'  )  J  dS  sin[7-7(s)] 

t'  rt' 


2iflhMsin(7 


wTTtJ'  dsJ 


dp  f (s)f (p)n(s)n(p) 


x  sin(7(s) -7]sin[7Cp) -7] 


(2.35) 


From  Eqs.  (2.31),  (2.32),  (2.34)  and  (2.35),  we  deduce  that  the  constants  b 


and  are  given  as 


-M7  Vx' 


b'  -  M7fjx 

o 


(2.36) 

(2.37) 


Also,  from  the  normalization  condition, 


do  - 


(2.38) 


From  Eqs.  (2 . 30) - (2 . 37)  and  (2.27'),  we  find  that 
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a(t,t' ) 

e(t.t') 

b(t.C') 

g(t.t') 

h(t.t') 

d(t.t') 


iM  ra  ^  1 


2¥  1 V 


+  7  COt(7-7' )  ] 


iM  .  n '  • 

—  [-^7  +  7'  coC(7-7') 


iM  _ 


XL. 


ft  sin(7-7') 

i  A 

ft  sin(7 


“I, 


ds  -  sin[7(s)*7] 

V7(s) 


*ii*ln<7-7-)  It.dS^7Slnl7'7<S)1 

-  •  ani,,  sinl7'7<s>l  Cdp 


x  sin[7(p)-7' 


(2.39) 

(2.40) 

(2.41) 

(2.42) 

(2.43) 

(2.44) 


Inserting  Eqs.  (2 . 39) - (2 .44)  in  Eq.  (2.11),  we  obtain  the  propagator  for  the  forced 
time-dependent  harmonic  oscillator  as 


K(x , t ;  x' t ' ) 


•  • 

M  (y'y) 


2mift  sin(7-7' ) 


.  iM  .n  2  tt'  ,  2 .  , 
x  exp[-  (f?  x  -  x'  )] 


f  i  M  r  •  2  •  2 

exp  \2ft  sin(7-7')  [(7X  +  7'x'  }  cos(y-7')  •  2Al‘ 


XX' 


zA 


„  X  [  ds  sin(7(s)-7'] 

Jt'  V7(s) 

+  2Ji~  x'  [  ds  sin[7-7(s)] 

Jt<  y;<S) 

-  'h  [  ds  zMMX  sin[7-7(s))  [  dp  sin[y(p) -7' ]]}  , 

JT  Jt'  JyCs)  Jt'  V7(P)  JJ 


(2.45) 


where  the  unpriraed  and  the  primed  variables  denote  the  quantities  which  are 
functions  of  time  t  and  t'  ,  respectively.  It  may  be  easily  verified  that  for  the 
case  where  w(t)  is  a  real  positive  constant  ,  we  have  rj  ( t)  -  1  and  7(t)  -  w^t, 
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and  the  propagator  of  Eq .  (2.45)  reduces  to  the  usual  expression  for  a  forced 
harmonic  oscillator . ^2 

III.  Wave  function 

We  now  rewrite  the  propagator  in  another  form  in  order  to  derive  the  wave 
function: 


K(x,t:  xt) 


m  Lthl 


2irity.  sin<7*7 


] 


x  exp 


“  {[?  *2  +  ^  x  I  dsJ^COSIl,'7<s)1. 


^7  x'2  +  jj/y'  x'  [  ds  -  cos  [7'  -  7(s)  ]iyi 

■"  J  y4<s>  JJJ 


x  exp 


^  cot(7-7*  ){[y-y  x  '  M  [  ds  7^  sin[7- 
■  J  J 7(s) 


7(s)  ] 


Jl'  X'  -  ^  f  ds  -Ui±  sin[7' -7(s)  ]12| 

M  J  y;<S)  J  j 


“  sinC^V,  {[^  X  •  M  I  dS  ^7  Sin[7'7< 

Jr*'  -  J  j 


s)  ] 


ds  sin[7' -7(s)  ]|YJ 

y-r(s)  JjJ 


X  exp[^5  -|cot(7-7')  U  ds  sin(7-7(s)] 


+  cot(7 


-'>[f 


dp  -iteX  sin[7'-7(p)]  2 


y-r(p) 


sin(7 


—  f 

-7')  J 


f' 


fl£l  . 


ds  sin(7-7^)  |  dp  ■— f -  sin[7’-7(p)] 

y4(P) 


1^777  {  ds-^sin[7-7(s)]  J  dp 

irU7  1  )  Jt,  Jy(s)  Jt'  J-, 


sin(l-T)  Jt, 

x  sin(7(p)-7' ]|j 


y  7(P) 


(3.1) 


Q 


M  U  .  .  J*  *  i(7*7f ) 

<:?  <i->  <*>  — 

1  -  e 


*  eXPS{[?  ^  +  M  A  X  j  dS  '“(T-TC*)]] 
x'2  +  M  ^'x'  J  dS  cos[7' -7(8)  ]jj>j 

exp[j{[y;'  *  •  M  |  ds  ^  sin[7"r(s)]]2 
77'X'  -  J  J  ds  sin[7'-7(s)]]2|‘] 


v;7(s) 


x  exp 

Li 


X[l-e 


-M _ 

■ 2i (7- 7 '  ) 


{  ^  X  '  M  f  ds  "7^  sin[7-7(: 
LL  M  J  M s) 


-  J  f 


ds  sin[  (7-7(s)  ]1  2 

V7(s)  J 


‘  2  ^  x  •  M  [  ds  sin(7-7(s)  ]1  [yyx' 

L  J  77(s)  jl 


-irds-^sin[7,.7(s)]iii 

J  J-y(s) 


e-U(t)eid(t') 


(3.2) 


‘  ^(t)  "  2KM  {cot(7-7')[[  ds  sin(7-7(s)]  2 

^  1  U  77(s)  J 

^  I 

+  cot(7-7')  j  ds  Oil  sin[7'-7(s>]  2  +  2  ~~ T~  7^ - — 

U  y7(s)  J  sm(7-7') 

x  [  ds  =^=  sin[7-7(s)]  [  dp  sin(7'-7(p)] 

J  y7(s)  J  y-^(p) 

-  [  ds  sin[7-7(s)]  f  dp  -1=^  sin(7<p)-7'] 

y7(s)  Jt-  y-^p) 


(3.3) 


at  us  introduce  Mehler's  formula, 


10 


2  2 


yi  -  z" 


2  2 

-  (X  +Y 


>  5  «  (X)H  ( 

L  0n,  n  n 


(3.4) 


“  o/t  x  +  jjj  f  ds  sin(7-7(s)]l 

^  J  v7(s)  J 


-  4  U'x'  -  J  f  ds  sin[7'-7(s)]l 

L  J  77<s)  j 


-1(7-7' ) 


Substituting  Eqs .  (3.4)-(3.7)  in  Eq.  (3.2),  we  obtain 


(3.5) 


(3.6) 


(3.7) 


K(x,t;  x',t)  -  )  0*(x,t)  0  (x',t')  , 

/  n  n 


(3.8) 


,  .  .  .  1  ,M~ir.  H.  H 

Vn(x.C)  -  [— ~  (  )  ]  exp 

2  n!  ’* 


fiMfi  2  2  ,♦  f*  f(s)  .  .  ,.1\ 

XPt2MU  X  ‘  M  77  X  J  dS  cos(7-(s)]J| 


exp  [ji  x  -  J  |  ds  sin[7-7(s)]]2} 

Hn  [X  *  J  dS  sin[7-7(sj]]> 

^  ^  L  m/7  j  y7(s) 


i[«(t)-(n+J)7(t)] 


(3.9) 


Moreover,  we  may  write 


0  ( x ,  t )  -  exp{i[i?(t)  -  (n+~)7(t)  ]  ]  0  (x,t)  , 

n  i.  n 


where 


VX-C)  ■  exp{2«[“  *'  M  J j  dS  co.(TT(.)]]} 


<3.10) 
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x  exp^JVy  X  -  j-j  |  ds  sin[7-7(s)  ]  j  2| 


*  Hn  '  M  I  ds  ■ln(T-T(«)l}] 


(3.11) 


In  Eq.  (3.10),  Che  wave  function  V>n(x,t)  is  merely  a  unicary  transformation  of 

i  (x,t),  and  thus  6  (x,t)  satisfies  all  the  properties  associated  with  (x,t): 
n  n  n 


dx  Tl>*ib  -  <m|n>  -  dx  4>*4>  -  & 

J  n*n  1  J  n^n  m,n 


(3.12) 


The  expectation  value  of  a  given  operator  0  is 


|o|n>  -  dx  V*  -  dx  <j>*0<f> 
11  J  n  n  J  n  n 


(3.13) 


IV.  Energy  expectation  values 

For  the  forced  time -dependent  harmonic  oscillator  system,  both  the  Hamiltonian 
and  Lagrangian  have  the  units  of  energy  but  depend  on  time.  We  must  therefore  find 
a  time- invariant  energy  operator.  If  /9(t)  is  a  particular  solution  of  Eq.  (2.3), 
we  have 


d  ° 

— -  (x-/9)  +  us  (t)  (x-/3)  -  0 

dt 


(4.1) 


and  from  Eqs.  (2.26)  and  (2.27')  we  note  that 


V  +  w2( t)  n  -  n2/r?3 


(4.2) 


From  Eqs.  (4.1)  and  (4.2),  we  get  the  following  expression  for  the  energy: 


1  *2  •  ••  M  •  •  2  M2  x  •  B  2 

E  "  2M  +  (Mr?x-^p)  (firi-rip)  +  ^(Pn-nP)  +  2  n 


(4.3) 


Because  Eq .  (4.3)  is  time  invariant,  we  can  use  it  for  the  quantum  mechanical 
energy  operator, 


“  <2*f;  ♦  1, 


-op  2K  Sx2  2' 


+  +  M^x)  +  Mf?72/3x  +  |  f?72/32  +  ^(Pb-V0)2 


(4.4) 
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Equation  (3.11)  now  simplifies  to  the  expression 


.  .  f  S  1 4  iux“+iAx  -hS^(x-S)^  u  ...  ,. 

V*'1'  -  7  8  8  H  [KW 

v 2  n'.Jir1 


-  f— H 

'■2nnl/f  J 


h  Ax"+Bx  „  .  r  .  ... 

e  5 (x-0) ]  , 


(4.5) 


where 


5  -  (M7/JO 


(4.6) 


M  2. 


(4.7) 


A ( t )  ~  u  Jy  f  ds  f cos(7-7(s) 
H  J  Jy(  s) 


(4.8) 


£(t)  -  [  ds  sin(7-7(s)l 

mVt  J  7t(s) 


(4.9) 


A  -  Lft  -  6  /2 


(4.10) 


B  -  LA  +  06‘ 


(4.11) 


tiere,  ;3(t)  is  a  particular  solution  of  Eq.  (2.3). 

In  order  to  evaluate  the  energy  expectation  E  -  <ra|E  |n>,  we  perform  the 

m ,  n  op 

following  calculations: 


x|n>  -  “TyT  [7n+l  |  n+l>  +  7n|n-l>] 


(4.12) 


x2|n>  -  ~~  (7(n+2) (n+1) | n+2>  +  (2n+l)|n>  +  7n(n-l) |n-2>] 


(4.13) 


p | n>  -  ^  2 (n+1 ) | n>  +  J  B|n>  +  ^(S+7)  72n|n-l> 


(4.14) 


2  2 

p2  |  n>  -  -jrf2  {~  7(n+2)  (n+1)  | n+2>  +  272  j1  7H+I|n>  +  2(A+^)  (2n+l)  |n> 


+  272  +  6B)7n|n-l>  +  2(4r+2A+62)7n(n-l) |n>] 


(4.15) 
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xp|n>  -  ^  V(n+2)  (n+1)  j  n+2>  +  -7=-  Jn+1 1 n+l>  +  [^(2n+l)+n]  | n> 


+  7=“  /H|n- 1>  +  (!-j+l)v/n(n+rTjn-2>} 

<5 


(4.16) 


px  |  n>  -  xp  |  n>  +  j  n> 


(4.17) 


Substituting  Eqs.  (4.6)-(4.17)  into  Eq.  (4.4),  we  directly  obtain  the  energy 
expectation  values  as 


En,n’  En-  2(^><2n+1) 


«h(n+~)  . 


(4.18) 


This  energy  expectation  value  is  a  time- invariant  quantity. 

V.  Uncertainty  values 

The  uncertainty  product  defined  as 

2  2*2  2  k 

(AxAp)^  ^  -  ( [ (<m|x  |n>  -  <m|x|n>  )  (<m|x  |n>  -  <mjx|n>  ) ] 

i2i  i  .  2  *  ,  2 i  ,  .  2  u  H 

x  [(<ra|p  j  n>  -  <m|p|n>  )  (<m|p  |n>  -  <ra|p|n>  )]') 


(5.1) 


Inserting  Eqs.  (4. 12) - (4. 15)  into  Eq .  (5.1),  we  obtain 


(axAp)n  n  -  U  -  ^!)h  <n*> 

7  1 


(5.2) 


( AxAp )  n+2  n  -  (1  +  7 7(n+2)  (n+1)  X 
’  7  n 


(5.3) 


(AxAp)n  n+2  -  (1  +  7 yn(n-l)  h  . 

7  n 


(5.4) 


VI.  Coherent  states  of  the  time -dependent  harmonic  oscillator 

First,  we  construct  the  creation  operator  a^  and  destruction  operator  a.  For 

a  forced  time -dependent  harmonic  oscillator,  it  is  not  possible  to  construct  a  and 
*  + 

a  .  but  ve  can  construct  a  and  a1  for  the  time  dependent  harmonic  oscillator.  From 
Eqs.  (4.12)  and  (4.14),  we  obtain 
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af-  (g)4  [(1  +  *7)  x  -  1^71 

m  My 


(6.1) 


a  -  (S^  ((1  *  i-2")  x  +  i2-] 


(6.2) 


From  Eqs.  (6.1)  and  (6.2),  we  can  represent  (x,p)  in  terms  of  (a  ,a)  as 


x  -  (-^-)4  (a^+a) 


(6.3) 


-  (^)4  [(*7  +  i)  af  +  (*7  -  i)a]  . 


(6.4) 


Also  from  Eqs.  (6.1)  and  (6.2),  if  [x,p]  -  i#  we  see  that 

[ a ^ , a ]  -  1  (6.5) 

Conversely,  from  Eqs.  (6.3)  and  (6.4),  if  [a^.a]  -  1  we  note  that  [x,p]  -  ijrf. 

The  coherent  state  can  be  defined  by  the  eigenstate  of  the  nonhermition 
15 

operator  a, 


a  a>  *»  a  a> 


(6.6) 


Let  us  find  the  coordinate  representation  of  the  coherent  state.  From  Eqs.  (6.2) 
and  (6.3),  we  have 


(S)h  [d  -  x  +  S  ^71  <x'  la>  "  a<x'  lQ>  • 

^  'll  My 

We  solve  this  equation  and  change  the  variable  x’  into  x  for  convenience, 


(6.7) 


|o>  -  N  exp  [IS(-1  +  i2?)  x2  +  (^S>4  a  x) 


We  choose  the  constant  of  integration  N  such  that 


dx  <x  a>  -  1 


(6.8) 


(6.9) 


Then,  we  find  the  eigenvector  of  the  operator  a  in  the  coordinate  representation 

j  x>  as 
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„  I  ^  /Mi.k  rMa,  .  ,fl_N  2  ,  ,2Mi  4  1|  ,2  1  ; 

<x|q>  -  (“£)  exp[^(-l  +  1^7)  x  +  (“g*)  a  *  -  -J | as  |  -  ?a 


(6.10) 


'll 


Nexc,  we  show  that  a  coherent  state  is  a  minimum  uncertainty  state.  From  Eqs . 

(6.3),  (6.4)  and  (6.6)  and  their  adjoints,  we  evaluate  the  expectation  values  of  x, 

2  2 

p,  x  and  p  in  the  state  |a>: 


<q|x|c»>  -  (~^-)^  (a  +a) 
2M 7 


(6.11) 


<a|p|a>  -  (^)^  [0*7  +  i)a*  +  -  i)ot] 


m 


ni 


(6.12) 


i  2 1  H  i  +2  2  -+-  ■+■  i 

<a  |  x  |  a>  —  —  <a  |  a  +a  +aa  +a  a  |  a> 

2My 


w  *2  2  * 

-JL-  (a  +a  +2aa  +1) 


2M7 


(6.13) 


•  • 


<q|p2|q>  -  ^  {(^7  +  i)2  a*2  +  (^7  -  i)2a" 
'll 


'll 


+  K3")2  +  1  ]  ( 2aa*+l )  )  . 


ni 


The  uncertainty  value  is 


o  o  u 

AxAp  -  [(<q|x  |a>  -  <a | x | a>) (<a | p  |a>  -  <a|p|a>)] 


(6.14) 


-  H/2  [l+(a7)2]'1 
'll 


(6.15) 


which  is  the  minimum  value  allowed  by  Eq.  (5.2) 


VII.  Results  and  discussion 

In  the  previous  sections,  we  have  obtained  the  propagator,  wave  function, 
energy  expectation  values ,  uncertainty  values  and  coherent  state  for  a  quantum 
forced  time -dependent  harmonic  oscillator.  These  quantities  represent  the  solution 
of  the  classical  equation  of  motion  for  the  time -dependent  harmonic  oscillator.  If 
we  set  f(t)  equal  to  zero,  then  our  solution  is  correct  for  the  time -dependent 
harmonic  oscillator.  Setting  w(t)  -  gives  results  for  the  forced  harmonic 


16 


oscillator.  For  the  case  where  f(t)  -  0  and  w(t)  -  w  ,  our  results  are  those  of 

o 

the  simple  harmonic  oscillator. 

For  the  explicit  time-dependent  system,  we  need  to  consider  the  quantum 

mechanical  operator.  In  our  work,  the  Hamiltonian,  Lagrangian  and  mechanical 

energy  have  the  units  of  energy,  but  these  are  not  time  invariant.  Yet,  in  order 

to  solve  macroscopic  physical  problems,  we  use  time -invariant  operators.  For  this 

reason,  we  have  derived  the  energy  operator  from  the  classical  equation  of  motion 

and  used  it  to  calculate  energy  expectation  values.  Our  energy  operator  is  similar 

1  2 

to  the  Ermakov-Lewis  invariant  operator.  ’  Our  quantum  energy  expectation  vlaues 
are  time- independent  quantities,  and  our  uncertainty  values  are  consistent  with 
Heisenberg's  uncertainty  principle.  Yet,  our  uncertainty  values  are  time 
dependent,  in  contrast  to  those  of  time- independent  systems. 

Since  it  is  not  possible  to  construct  a  coherent  state  for  the  forced  time- 
dependent  harmonic  oscillator,  we  have  constructed  it  for  the  time  -  dependent 
harmonic  oscillator.  In  general,  the  coherent  state  is  a  minimum  uncertainty 
state,  which  is  also  true  for  our  system. 

Time -dependent  systems  are  observed  in  various  physical  experiments.  Two 
general  types  of  such  systems  are:  that  which  is  formed  through  its  own 

environmental  conditions,  and  that  which  is  formed  when  external  forces  are  added. 
In  regard  to  the  second  type,  various  experiments  are  being  carried  out  to  see  how 
an  applied,  time -dependent  electric,  magnetic  or  other  field  can  alter  the  physical 
properties  of  materials  such  as  semiconductors  and  superconductors.  Experiments 
show  that  a  system  becomes  time  dependent  when  a  time-dependent  electric  or 
magnetic  field  (such  as  a.c.)  is  applied.  However,  obtaining  the  quantum 
mechanical  solution  by  a  direct  method  is  not  easy  mathematically.  One  way  of 
obtaining  a  solution  is  to  use  the  propagator  method  as  indicated  in  this  paper, 
where  the  relevant  equations  are  those  of  a  time-dependent  harmonic  oscillator. 

Our  results,  which  are  exact  for  one  dimension,  can  be  extended  to  two  or  more 
dimensions,  and  they  can  also  be  applied  to  time -dependent  macroscopic  systems. 
One  example  of  an  extension  to  two  dimensions  would  be  to  solve  the  motion  of  a 
quantum  electron  in  a  time -dependent  magnetic  field. 

Acknowledgments 

KWY  and  CIU  acknowledge  financial  support  from  the  Korea  Science  and 
Engineering  Foundation  and  the  Basic  Science  Research  Institute  (BSRI),  Ministry  of 
Education  1989,  Republic  of  Korea.  This  research  was  also  supported  by  the  Office 
of  Naval  Research  and  the  National  Science  Foundation  under  Grant  CHE- 9016789. 


17 


References 

1.  H.  R.  Lewis,  Jr.,  Phys.  Rev.  Lett.  18.,  510,  636  (1967). 

2.  H.  R.  Lewis,  Jr.,  J.  Math.  Phys.  9,  1976  (1968). 

3.  P.  G.  L.  Leach,  J.  Math.  Phys.  10,  1902  (1977). 

4.  J.  R.  Burgan,  M.  R.  Felx,  E.  Fijalkow  and  A.  Munier,  Phys.  Lett 
(1979)  . 

5.  J.  G.  Hartly  and  J.  R.  Ray,  Phys.  Rev.  A  24,  2873  (1981). 

6.  P.  Camiz,  A.  Gerardi,  C.  Marchioro,  E.  Presutti  and  E.  Scacciatelli , 
Phys.  12,  2040  (1971). 

7.  D.  C.  Khandekar  and  S.  V.  Lawande,  J.  Math.  Phys.  16.,  384  (1975). 

8.  A.  D.  J annus sis,  C.  N.  Brodimas  and  A.  Streclas,  Phys.  Lett.  A  74,  6 

9.  C.  I.  Um,  K.  H.  Yeon  and  W.  H.  Kahng,  J.  Phys.  A  20,  611  (1987). 

10.  K.  H.  Yeon,  C.  I.  Um  and  T.  F.  George,  Phys.  Rev.  A  36,  5287  (1987). 

11.  K.  H.  Yeon,  C.  I.  Um,  W.  H.  Kahng  and  T.  F.  George,  Phys.  Rev.  A 

(1988)  . 

12.  R.  P.  Feynman  and  A.  R.  Hibbs ,  Quantum  Mechanics  and  Path  Integrals 
Hill,  New  York,  1965). 

13.  G.  J.  Papadopoulous ,  Phys.  Rev.  D  H,  2870  (1975). 

14.  A.  Erdelyi,  Higher  Transcendental  Functions,  Vol.  2  (McGraw-Hill, 
1953),  p.  194. 

15.  W.  H.  Louisell,  Quantum  Statistical  Properties  of  Radiation  (Wiley, 
1973)  . 


.  A  74,  11 


J.  Math. 


(1979). 


1£.  6224 


(McGraw 


New  York, 

New  York, 


FY91  Abstracts  Distribution  List  for  Solid  State  f <  Surface  Chemistry 


Dr.  Paul  Ballentine 
CVC  Products 

525  Lee  Rd. ,  P.  0.  Box  1886 
Rochester,  NY  14627 


Dr.  Andrew  Baronavski 
Code  6111/Chemistry  Division 
Naval  Research  Laboratory 
Washington,  D.C.  20375-5000 


Dr.  John  Crowell 
Department  of  Chemistry 
University  of  California 
La  Jolla,  CA  92093 


Dr.  Frank  DiSalvo 
Department  of  Chemistry 
Cornell  University 
Ithaca,  NY  14853 


Dr.  John  Eyler 
Department  of  Chemistry 
University  of  Florida 
Gainesville,  FL  32611 


Dr.  Paul  Barbara 
Department  of  Chemistry 
University  of  Minnesota 
Minneapolis,  MN  55455-0431 


Dr.  R.P.H.  Chang 
Dept.  Mat.  Sci.  &  Eng. 
Northwestern  University 


Dr.  Mark  D' Evelyn 
Department  of  Chemistry 
Rice  Univ.,  P.  0.  Box  1892 
Houston,  TX  77251 


Dr.  Arthur  Ellis 
Department  of  Chemistry 
University  of  Wisconsin 
Madison,  WI  53706 


Dr.  James  Garvey 
Department  of  Chemistry 
SUNY/Buffalo 
Buffalo,  NY  14214 


Dr.  Robert  Gedridge 
Code  3854/Chemistry  Division 
Naval  Weapons  Center 
China  Lake,  CA  93555-6001 


Dr.  Paul  Hansma 
Department  of  Physics 
University  of  California 
Santa  Barbara,  CA  93106 


Dr.  John  Hemminger 
Department  of  Chemistry 
University  of  California 
Irvine,  CA  92717 


Dr.  Leonard  Interrante 
Department  of  Chemistry 
Rensselaer  Polytech.  Inst. 
Troy,  NY  12181 


Dr.  Steven  George 
Department  of  Chemistry 
Stanford  University 
Stanford,  CA  94305 


Dr.  Robert  Hamers 
Department  of  Chemistry 
University  of  Wisconsin 
Madison,  WI  53706 


Dr.  Charles  Harris 
Department  of  Chemistry 
University  of  California 
Berkeley,  CA  94720 


Dr.  Roald  Hoffmann 
Department  of  Chemistry 
Cornell  University 
Ithaca,  NY  14853 


Dr.  Eugene  Irene 
Department  of  Chemistry 
University  of  North  Carolina 
Chapel  Hill,  NC  27514 


Dr.  Zakya  Kafafi 
Code  6551 


Naval  Research  Laboratory 
Washington,  DC  20375-5000 


Dr .  Fred  King 

Dept,  of  Chem,  PO  Box  6045 
West  Virginia  University 
Morgantown,  WV  26506-6045 


Dr.  Stephen  Lieberman 
Marine  Evn.  Branch/Code  522 
Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 


Dr.  Horia  Metiu 
Department  of  Chemistry 
University  of  California 
Santa  Barbara,  CA  93106 


Dr.  Daniel  Neumsrk 
Department  of  Chemistry 
University  of  California 
Berkeley,  CA  94720 


Dr.  Larry  Kesmodel 
Department  of  Physics 
Indiana  University 
Bloomington,  IN  47403 


Dr.  Max  Lagally 
Dept.  Metal.  S.  Min.  Eng. 
University  of  Wisconsin 
Madison,  WI  53706 


Dr.  M.C.  Lin 
Department  of  Chemistry 
Emory  University 
Atlanta,  GA  30322 


Dr.  Larry  Miller 
Department  of  Chemistry 
University  of  Minnesota 
Minneapolis,  MN  55455-0431 

Dr.  David  Ramaker 
Department  of  Chemistry 
George  Washington  University 
Washington,  D.C.  20052 


Dr.  Gary  Rubloff 

IBM  T.J.  Watson  Res.  Center 

P.  0.  Box  218 

Yorktovm  Heights,  NY  10598 


Dr.  Richard  Smalley 
Department  of  Chemistry 
Rice  Univ.,  P.0.  Box  1892 
Houston,  TX  77251 

Dr.  Galen  Stucky 
Department  of  Chemistry 
University  of  California 
Santa  Barbara,  CA  93106 


Dr.  William  Unertl 

Lab.  for  Surface  Sci.  &  Tech. 

University  of  Maine 

Orono,  ME  04469 


Dr.  John  Weaver 

Dept,  of  Chem.  &  Mat.  Sciences 
University  of  Minnesota 
Minneapolis,  MN  55455 

Dr.  Howard  Schmidt 
Schmidt  Instruments,  Inc. 

2476  B^lsover,  Suite  234 
Houston,  TX  77005 


Dr.  Gerald  Stringfellow 
Dept,  of  Materials  Sci.  &  Eng. 
University  of  Utah 
Salt  Lake  City,  UT  84112 


Dr.  R.  Stanley  Williams 
Department  of  Chemistry 
University  of  California 
Los  Angeles,  CA  90024 

Dr.  Aaron  Wold 
Department  of  Chemistry 
Brown  University 
Providence,  RI  02912 


Dr.  John  Yates 
Department  of  Chemistry 
University  of  Pittsburgh 
Pittsburgh,  PA  15260 


Dr.  Robert  Whetten 
Department  of  Chemistry 
University  of  California 
Los  Angeles,  CA  90024 


Dr.  Nicholas  Winograd 
Department  of  Chemistry 
Pennsylvania  State  University 
University  Park,  PA  16802 


Dr.  Vicki  Wysocki 
Department  of  Chemistry 
VA  Commonwealth  University 
Richmond,  VA  23298-2006 


Dr.  H.  Tachikawa 
Department  of  Chemistry 
Jackson  State  University 
Jackson,  MI  39217-0510 


i  Dr.  Terrell  Vanderah 
Opt.  Elec.  Matls/Code  3854 
Naval  Weapons  Center 
China  Lake,  CA  93555 


Dr.  Brad  Weiner 
Department  of  Chemistry 
University  of  Puerto  Rico 
Rio  Piedras,  PR  00931 


Dr.  Paul  Weiss 
Department  of  Chemistry 
Pennsylvania  State  University 
University  Park,  PA  16802 


